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Abstract 

Quantum mechanical behavior of coupled A^-kicked rotators is studied. In the 
large N limit each rotator evolves under influence of the mean-field generated by 
surrounding rotators. It is found that the system spontaneously generates classical 
chaos in the large limit when the system parameter exceeds a critical value. 
Numerical simulation of a quantum rotator coupled to a classical rotator supports 
this idea. 
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1 Introduction 



In spite of numerous efforts to define quantum chaos, there seems to exist many counter 
examples of genuine chaos in quantum systems. Among those the quasi-periodic time 
evolution of isolated bound systems and the limited diffusion due to the momentum 
localization in the periodically kicked quantum rotator are typical examples Al- 
though they are rather persuasive to give up searching quantum chaos, we still reserve a 
good reason to believe the seeds of chaos should exist in quantum system which grow to 
classical chaos. One way to observe it is in the time evolution of the Wigner function; 

When the quantum Liouville operator L is expanded in terms of the Planck constant as 

L = Lo + n^Li + ■ ■ • , (2) 
the lowest order term generates the classical evolution, 

W] = Y.{hjY% (3) 

i 

where {if, /}'•*'' represents the Poisson bracket with respect to the i-th canonical variables 
{qi,Pi) of a system with the Hamiltonian H{{qi}, {pi})- Therefore, it is natural to search 
the origin of chaos in the region of ~ 0, i.e. the semi-classical limit. In fact, recent 
intensive investigations [^,0 have revealed rich structures in this limit. 

There is another way to observe the emergence of chaos. If it exists in the system 
which involves many degrees of freedom. Since the most characteristic property of a 
macroscopic system is the large number of degrees of freedom involved in describing the 
motion, it is natural to consider a system with infinite number of degrees of freedom as the 
realistic classical limit. In this case we expect that quantum mechanical phases of each 
degree of freedom interfere mutually, and are wiped out in chaotic environment resulting 
the classical behavior of the Wigner function, i.e. the classical distribution function 0. 

When the number of degrees of freedom involved in a system tends to infinity, the 
mean-field approximation often becomes exact due to the large number theorem. In such 
a case the Wigner function can be approximated in a product form O, 

N 

f{{qi}AP^hi)^]lf^iQ^,P^■^i), (4) 
i=l 
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of functions corresponding to the sub-systems. The Wigner function for the i-th sub- 
system then evolves as 

^ = Lf\U, (5) 

with the appropriately defined mean- field operator, Lf\ In general, this operator depends 
on the temporal distributions of surrounding sub-systems. Thus the mean-field Liouville 
equation (^) becomes non-linear and time- dependent. It is our expectation that the non- 
linearity and the time- dependence induce chaos, and chaos destroys quantum phases. 

In the next section we shall describe time evolution of the Wigner function for a 
system of coupled kicked rotators, and the mean-field approximation is introduced. In 
the third section the mean-field equation is replaced by an equivalent system of a quantum 
kicked rotator coupled with a classical kicked rotator representing the environment. And 
numerical simulation of this system is examined. The idea of spontaneous decoherence is 
discussed based on the numerical calculation in the last section. 

2 Large N limit and the statistical mean-field ap- 
proximation 

Let us examine the scenario more explicitly by taking the coupled iV-kicked rotator model 
1^ as an example. The Hamiltonian is assumed to be 

N ^ N 
i=l ^ ^ ij^j 

where 

oo 
n=—oo 

is the periodical kick. In this model couplings of kicked rotators enter through off- 
diagonal kinetic terms. Couplings through the kick term [^] may be converted to this 
form through a linear transformation of angle variables. The classical equation of motion 
is known as the standard map which is the mapping between canonical coordinates just 
after kicks at times t and t -|- 1, 

jm) _ „W , x'Ji+i) {t: integer) (7) 
Pi - Pi + K sm Oi 
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We shall choose the magnitudes of matrix elements of the kinetic term as 

'« ~ °(^)- (8) 
~ 0{^) for i # j, 

where the exponent u will be determined later so that the system is non-trivial and 
converging in the large limit. 

Evolution of the Wigner function is obtained by solving the Schrodinger equation. 
The relation between Wigner functions just after kicks at times t and i + 1 is called the 
quantum mapping, and it can be written as 

f'-''\{o^}Ap^}) = uvf'\{e,},{p,}) (9) 

with the mapping operators U and V defined by 

C/ = exp(-^AiSin^,-^), (10) 

and 

V^eM-EtijPi^)- (11) 



Here, the operator appeared in U, 



D 1 , h_a_ _h8_^ 

_ = _(e2ap _e 2ap), (12) 

is a difference operator, and quantum effects in the evolution equation enter only through 
this operator. Since it approaches to the differential operator in the classical hmit as 

§-^Hp) = > + ^)-/(P-^)l 

~ ^/(f),(R-o) (13) 

the quantum mapping reduces to the well-known standard map of classical distribution 
function, 

f'^-'Hrn, {P^}) ^ f^'\{e^ -J2t^JiPJ ' A, siu^,)}, {p, - A.sin^,}) (14) 

j 

in the — > limit. 
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Although the difference between classical and quantum evolutions looks very small 
from eqs.(|TO[) and([T3|) at a glance, the consequence of actions of the operator U on the 
exponential function may be instructive. In quantum mechanical case it acts as 

with 

while in the classical limit it reduces to the standard shift operator along the momentum 
axis, 

g-'^l^gifcp = gifc(p-a)^ (17) 



i.e. the operator behaves as the usual shift operator only for small wave number 
components. For large wave numbers shifts are suppressed by the factor 1/k and direction 
of the shift alternates. This implies that rapid variations of the Wigner function along 
the p-axis are suppressed, and consequently chaotic folding of the distribution function 
terminates at fine scale of the order h. 

Now, let us consider the large N limit where we expect that the mean-field approxi- 
mation works well. Let us suppose at a certain time t the Wigner function is written in 
a product form, 

N 

f^'\{0^},{P^}) = I[f^'\O^.P^)■ (18) 
i=l 

Integrating the evolution equation(^ over all phase volume except the i-th canonical 
coordinates, we obtain the equation for the reduced Wigner function of the i-th sub- 
system: 

D d 
ft^^\Oi,Pi) = exp(-AiSin6'i— — )exp(-tiiPi— ) 

Dpi d&i 

d 

Mi) 

where the reduced Wigner function has been defined at t + 1 as 



X 



I n rf^,rfp,exp(- (19) 



^t'\Q^.P^) = / n dQ,dv,f'^'^'\{Q,}APi\). (20) 

which will turn out to be the solution of the equation which we may call the statistical 
mean-field equation evolved from fi^\Oi,pi) of eg . (|l^) . 
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Expanding the exponential function in the integrand of the r.h.s. of eq.(]T9|) we obtain 

d 



+ \ E uMp.)f{p,)f^,+Wtl{p])f^, 



for the first and second order terms, where we have written 

{p^)f = Jde,dp,p]ff{e„p,). (21) 

Writing back this expression to the exponential form we have 



exp(-E^.(p.)f^+E4-f^)i^'^ 



for the first two terms of the cumulant expansion. Here, cr^-*^ is the momentum variance 
of the j-th sub-system, 

4' = (j^)?-m?r- (22) 

We have chosen that the coupling strength tij{i ^ j) to be 0{N~'^), and thus the n-th 
cumulant {i.e. the term involving ^) is of the order of N^^™ . The parameter u classifies 
nature of the system into the following cases in the N ^ oo limit: 

i) z/ > 1, any cumulant converges to and the system becomes a set of independent 
kicked rotators. 

ii) = 1, only the first cumulant tij{Pj)f^ survives and the ordinary mean-field 
Hamiltonian, 

hi = hi + Y^ PiUj {pj)f , (23) 

drive the reduced Wigner function //*■*. 

iii) I < z/ < 1, the ordinary mean-field theory is valid assuming J2j(j^'i)tij{Pj)f^ re- 
mains finite. By transforming the system in a rotating frame with an appropriate 
momentum p, the fluctuation tij{{Pj)f^ — p) can be kept finite. 

iv) v = \, the statistical mean-field theory is valid providing J2j{j^i)tij{{pj)f^ — p) kept 
finite. 
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v) z/ < |, the second cumulant diverges and the system becomes unphysicaL In this 
case it does not make any sense to consider the reduced Wigner function. 

The statistical mean-field equation is given by 

D d 
fi^^'^\di,Pt) = exp{-XiSmei——)exp{-tiiPi—) 

Dpi dOi 

exp(- E i^^^P^)f^. + E (24) 

Since evolution operator of the i-th subsystem involves only the i-th canonical coordinates 
product form of the Wigner function will remains at t + 1: 

N 

f'^'\{e^},{V^}) = \{ft'\e^.V.). (25) 

i=l 

The proof that the product eq. (|25| ) is indeed the solution of the evolution equation (|) is 
given in ref. 0. 

For studying asymptotic behavior of the momentum distribution it is more convenient 
to consider the Fourier transform of the reduced Wigner function in the statistical mean- 
field approximation, 

where mi and fcj are the Fourier parameters corresponding to the variables 6i and pj, 
respectively. The momentum distribution function for the i-th sub-system is then readily 
expressed as 

ff\p)= r ^e-^''/^jT{0,k). (26) 
The evolution equation in this representation is written as 

f~r\m., k.) = E ^.-k(x ^)e-^'''-^-^''^'f^i< - (27) 
where Jm{k) is the Bessel function of the order m, and we have defined 

gf^ = E t^Ap^)f\ (28) 

= E • (29) 
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For the asymptotic behavior of sf'*, there will exist at least two cases depending on 
the nature of environment: 

crf ■* —>■ const, (regular) 

dit. (chaotic) (30) 



When the environment is chaotic the sum over m, in eq.(27) will be dominated asymp- 
totically by the mj = components. The asymptotic form of the evolution equation then 
becomes 

ft'\m.M - Jn.A'^-^sin^)f7\0,k,), (31) 
whose solution is now trivial. Especially the = components are given as 

ff^{0,k^ ~ {M'^-^sm^)yfT\0,h). (32) 



Inserting this solution into eq. ( P6|) together with help of the cental limit theorem we obtain 
the asymptotic form of the momentum distribution function as 

#^(p)~[^]^e"%, (33) 

which is the same function as the classical distribution of the chaotic limit. The m = 
dominance suggests that when the environment become chaotic the system starts behaving 
classically due to the phase decoherence. 



3 An equivalent model and numerical results 

We have seen that the A^-coupled kicked rotator will fall into chaos in the large limit 
in cases when the statistical mean-field approximation works well. Under those circum- 
stances it becomes possible to study, instead of treating whole system at once, just by 
picking up any one of the rotators as an object system which we simply call the system(S') 
under influence of other rotators as the environment (i?). Interaction between 5* and E 
enters as an external field with strength proportional to the state averaged momentum in 
our model. Although the mean-field approximation makes it possible to perform the large 
calculation as many as ~ 100 — 1000 by using large scale parallel processors, we leave 
such huge simulations as the future project. Instead, we assume the mean-field strength 
(PSP is replaced by a classical system with appropriate initial ensemble. We choose the 
ensemble so that it become equivalent to the statistical mean-field theory. 
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Now, let us consider a quantum kicked rotator(S') with the Hamiltonian, 

Hs = ^p' + A{t)Xs cos e + g^'^'p (34) 
under the influence of classical environment through g^'^\ This system evolves as 

ft'\e,p) = exp(-A5sin^-^)exp(-p^)exp(-^7(*)^)/W(e,p), (35) 
or by the Fourier transform it reads 

fT'\m, k)=Y: Jm-m'{^^ Sin ^)e~-'<^'V^J\m\ k - m'). (36) 

m' 

After taking average over the initial ensemble the evolution equation is written as 

E[ft'\m. k)] = E Jm-m\^^ sin ^)E[e-'^'^''' fl\m\ k - m% (37) 

m' 

where we write the ensemble averaged quantities by £'[ ]. We assume the ensemble of 
mean-field couplings {g^^^*} to be Gaussian with the correlation, 

i^[^W^(*')]_i5[^W]i5[^(t')] = 5^^,^W. (38) 

In this case the correlation between g^*^ and /'•*^ vanishes, and we have 

E[ft'\m, k)] = J„,(^ sin |A:)e— '^[^'"l-^-'^-"'E[/^*)(m', k - m')]. (39) 

m' 

This equation is equivalent to the statistical mean-field equation (pT]). 

Our next task is to give g^^^ explicitly in order to solve the evolution eq.(^3]). One of the 
possible choice of g^^^ which possesses both properties eq.(pO|) and (|38|) is the momentum 
generated by the classical kicked rotator, 

g{t+i) = (;W + A£;sin0(*+i) 



with a set of initial values g^^^) as the ensemble. Iterating eq. (|35D together with 
eq.(^OD the behavior of S will show the quantum-classical transition depending on the 
environment. 

Before examining results of the numerical calculation let us remark the consistency of 
the choice of g^^\ The momentum distribution of E in the mean-field approximation is 
given by 

f^\pE) = / n dp^PE- E p.) n /f fe)- (41) 
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If each rotator is in a chaotic state its momentum distribution is known to be 

ff\p)(xe~^ (42) 

asymptotically. In this case fE\pE) also has the same form with the diffusion constant 
being 

dE=Y. tyj, (43) 

which is C(l) for our choice of tsj ~ 

We shall now examine numerical results which support the scenario described in the 
last section. We solve the Schrodinger equation with the Hamiltonian (p^ with g^^'> 



generated from the standard map (^Ol). At each time step we measure the momentum 
fluctuation of S, 

averaged over the initial ensemble which consists typically 20 members for each set of 
parameters A5, A^; and fi. As in the classical rotator chaos is measured by the diffusion 
constant defined by the rate of momentum fluctuation. 

In order to see the effect of chaotic environment the system is maintained in chaos 
by choosing = 6.0, and varying the environment A^; (Fig.l). The coupling strength is 
chosen to be yU = 0.01. As we have expected the energy diffusion starts precisely when 
the environment become chaotic, i.e. Xe ~ 1.0 . To make sure that the increase of kinetic 
energy is not due to the direct penetration of energy from E via coupling, we made the 
similar calculation by keeping the system in the regular regime, A5 = 0.8. No energy 
diffusion has been observed in this case. 

Choice of the coupling strength fi = 0.01 comes from the requirement that the property 
of S should not be altered significantly by the environment. Weakness of the coupling 
effect is clear if we compare diffusion constants of an isolated classical rotator which is 
18.0 for A = 6.0. No significant difference is observed when the coupling strength is less 
than 0.01. Although the effect of environment is small on the classical system, it has 
the essential action on the quantum system. It is seen in Fig.2 where S and E are kept 
in chaos (Ag = A^; = 6.0) and the coupling strength is varied. Above the lower critical 
coupling /i ^ 0.001 the environment effect appears in the diffusion constant of S. 

Next, we observe the transition of S from the regular motion to chaos under the 
influence of the chaotic environment with keeping A^; = 6.0 (Fig.3). At the classical 
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critical strength A = 1.0 the quantum system persists to be diffusion- less. It is only after 
A ~ 3.0 the system shows diffusion. Change of the evolution from quantum to classical 
by the effect of E can also be seen in the momentum distribution. Figs. 4 (a) and (b) 
correspond to the regular motion and the chaotic regime, respectively. In the regular case 
the momentum localization persists during the simulation time up to t < 300, while in 
the chaotic case the momentum distribution changes following the classical diffusion law 
described by eq. (|53D . 

4 Spontaneous decoherence -A conclusion- 
it is well-known that a single quantum kicked rotator shows only limited energy diffusion. 
The limited diffusion is considered to be generic in quantized systems with a few degrees of 
freedom whose motion show chaos in the classical evolution. We have searched the source 
of classical diffusion in systems with infinite number of degrees of freedom. These system 
can be treated in the statistical mean-field approximation, where we can pick any one of 
the sub-system as an object system and the rest as the environment. Energy diffusion 
has been observed in the quantum system when the environment is in classical chaos. 
Momentum distribution of the quantum system is also same as one which we normally 
observe in a classical chaotic system. Since each building block of the environment is the 
same type of quantum system, the transition from the quantum behavior to the classical 
behavior seen in the diffusion should occur spontaneously. 

One may wonder then where is the seed of chaos in linear dynamical system like 
quantum mechanical system. It is certainly very difficult to point out the cause of chaos 
unlike classical dynamical system where the existence of homoclinic point means chaos. 
Although we have no rigorous proof for the existence of chaos in quantum systems in 
the limit of infinite number of degrees of freedom, the exactness of the certain class of 
mean-field approximation in the large N limit suggests that the non-linearity is secretly 
inherited in quantum system with infinite number of degrees of freedom. It is the non- 
linear mean-field equation which possesses the source of chaos within itself. 
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Figure 1: Diffusion constants of S fixing the kick strength in the chaotic regime(A5 = 6.0) 
and varying to change the environment from regular to chaos. 
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Figure 2: Diffusion constants of S fixing the kick strength of S and E both in the chaotic 
regime(A5 = A^ = 6.0), and varying the coupling 
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Figure 3: Diffusion constants of S fixing the kick strength of E in the chaotic 
regime(A£; = 6.0), and varying A5. 
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Figure 4: Time evolutions of the momentum distributions of S setting the kick strength 

of E in the chaotic rcgimc(As = 6.0) and keeping the system in (a) the regular motion 
{\s = 0.8), or (6) the chaotic motion (As = 6.0). Both of (a) and (6) are given by 100 
ensembles average. 
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